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Two-Level Rectilinear Steiner Trees 


Stephan Held* Nicolas Kammerling* 


Abstract 

Given a set P of terminals in the plane and a partition of P into k subsets Pi,..., Pk, 
a two-level rectilinear Steiner tree consists of a rectilinear Steiner tree Tj connecting 
the terminals in each set Pi {i = 1,... ,k) and a top-level tree Ttop connecting the trees 
Ti,... ,Tk. The goal is to minimize the total length of all trees. This problem arises 
naturally in the design of low-power physical implementations of parity functions on a 
computer chip. 

For bounded k we present a polynomial time approximation scheme (PTAS) that is 
based on Arora’s PTAS for rectilinear Steiner trees after lifting each partition into an 
extra dimension. 

For the general case we propose an algorithm that predetermines a connection point 
for each Tj and Ttop (* = 1,... , A:). Then, we apply any approximation algorithm for 
minimum rectilinear Steiner trees in the plane to compute each Tj and Ttop indepen¬ 
dently. 

This gives us a 2.37-factor approximation with a running time of 0(|P| log |P|) suit¬ 
able for fast practical computations. The approximation factor reduces to 1.63 by 
applying Arora’s approximation scheme in the plane. 


1 Introduction 

We consider the two-level rectilinear Steiner tree problem (R2STP) that arises from an application 
in VLSI design. Consider the computation of a parity function of k input bits using 2-input XOR- 
gates. Due to the symmetry, associativity, and commutativity of the XOR function, this can be 
realized by an arbitrary binary tree with k leaves, rooted at the output, by inserting an XOR- 
gate at every internal vertex [13]. Throughout this paper we consider the parity function as a 
placeholder for any fan-in function of the type xi 0 x 2 o • • • ox^, where o is a symmetric, associative, 
and commutative 2-input operator, i. e. o G {©,V, A}. 

On a chip such a tree has to be embedded into the plane and all connections must be realized 
by rectilinear segments. If each input and the output are single points on the chip, a realization of 
minimum length and thus power consumption is given by a minimum length rectilinear Steiner tree. 
This is a tree connecting the inputs and the output by horizontal and vertical line segments using 
additional so-called Steiner vertices to achieve a shorter length than a minimum spanning tree. At 
each Steiner vertex of degree three an XOR-gate is placed. Higher degree vertices can be dissolved 
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Figure 1: On the left, we have two inputs pi and p 2 and a single output p^. The XOR-gate should 
be placed at the median of the three terminals. If the inputs have the side outputs p'l 
and P 2 , the XOR-gate should be placed at ps, saving the horizontal length. 


into degree three vertices sharing their position. Figured] shows an example of an embedded parity 
function on the left. 

In practice input signals may be needed for other computations on the chip and thus delivered 
to other side outputs. Similarly, the result may have to be delivered to multiple output terminals. 
Thus, each input and its successors and the output terminals must be connected by separate Steiner 
trees as well. These trees are then connected by a top-level Steiner tree into which the XOR-gates 
will be inserted. Considering the additional terminals allows to construct a potentially shorter 
top-level and two-level Steiner tree as shown in Figured) on the right. Algorithms ignoring the side 
outputs cannot guarantee an approximation factor better than two, as we will see in Section [2] 

This motivates the definition of the minimum two-level rectilinear Steiner tree problem, where 
we are given a set P C of n terminals and a partition of P into k subsets Pi,..., P^- 

A two-level rectilinear Steiner tree T = {Ttop, Ti,..., Tfe) consists of a Steiner tree Ti for each i € 
{1,..., k} connecting the terminals in Pi and a (group) Steiner tree Ttop connecting the embedded 
trees {Ti,... ,Tk}. We call Ttop the top-level tree. Note that all trees are allowed to cross. The 
objective is to minimize the total length of all trees 

k 

1{T) ■.= Y,KTi)+l{Ttop), 

i=l 

where 1{T') := II® “ ^lli .^i-length of a Steiner tree T'. 

For each i G {1,..., fe} the top-level tree and Tj intersect in at least one point. We can select one 
such point qt G Ttop H Ti and call it connection point for Ti and Ttop- Then Ttop is a Steiner tree for 
the terminals {gi,..., qk} and each Ti is a Steiner tree for Pi U {qi}. 

Obviously, this problem is JVP-hard as it contains the minimum rectilinear Steiner tree problem 
in two ways: if = 1 or if |Pj| = 1 for i G {1,..., k}. The problem is JVP-complete, because there 
is always an optimum solution in the Hanan grid [9] of P. This simple fact will arise later as a 
side-result in Corollaries EH and EH 

Designing the top-level tree as a stand-alone problem is hard. If all subtrees Ti {i G {1,..., k}) 
are hxed, Ttop cannot be approximated to arbitrary quality, as the group Steiner tree problem 
for connected groups in the Euclidean plane cannot be approximated within a factor of (2 — e) 
dH. However we are in a more lucky situation as we can tradeoff the lengths of bottom-level and 
top-level trees. 

To the best of our knowledge the two-level rectilinear Steiner tree problem has not been considered 
before despite its practical importance [mils]. It is loosely related to the hierarchical network 
design problems mum or multi-level facility location problems mu- However, those problems 
are structurally different, typically considering problems in graphs, and do not apply to our case. 
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In [T3], ordinary rectilinear Steiner trees were used to build power efficient fan-in trees, when 
each input and the output consists of a single terminal. In practice designers are also interested in 
the depth of the constructed circuit m- However, for finding good power versus depth tradeoffs a 
better understanding of short solutions is an essential prerequisite and the aim of our work. 

1.1 Our Contribution 

In Section[2]we show even for arbitary metrics that the naive approach of picking a random terminal 
from each partition as a connection point to the top-level tree and building the bottom-level trees 
and top-level as separate instances gives a 2a-factor approximation, where a is the approximation 
factor of the used minimum Steiner tree algorithm. 

Then in Section [3] and 0] we focus on rectilinear instances. In Section [3] we show how to lift our 
instance into an equivalent (2 -|- A:)-dimensional rectilinear Steiner tree instance. If the number k of 
partitions is bounded by a constant, we obtain a PTAS by applying Arora’s PTAS for rectilinear 
Steiner trees [2]. 

As our main result we improve the approximation guarantee for unbounded k from (2-|-e) to 1.63 
in Section 01 Using spanning tree heuristics this approach turns also into a fast practical algorithm 
with running time O(nlogn) and approximation factor 2.37. 

2 Simple Bottom-Up Construction 

A simple bottom-up approach, which works for any metric space, is to compute a Steiner tree Tj 
for Pj (i = 1,..., fe). In each Tj we fix a connection point qi € Pi arbitrarly, compute a Steiner tree 
Ttop for {qi ,..., qk}, and return T = {Ttop, Ti, • • •, Tfc). 

Theorem 2.1. The simple hottom-up approach is a 2a-factor approximation algorithm for the 
minimum two-level Steiner tree problem, if we use an a-factor approximation algorithm for the 
minimum Steiner tree problem as a subroutine. 

Proof. Let T be the two level Steiner tree computed by the simple bottom-up approach and let 
T* = (Tj*p, Tf,... ,T^) be a minimum two-level Steiner tree. Let be q* G n Tf the connection 
point of the optimum two-level Steiner tree. Since T* is a Steiner tree on {g*} Li Pi, we have 
dist{q*,qi) < 1{T*). Thus, 

k k k 

1{T) =l{Ttop)+'^l{Ti) < a-l{Tf^)+aJ2dist{q*,qi) + '^a-l{T*) 
i=l i=l 2=1 

k 

< a ■ + 2a Y, 1{T*) < 2a-l{T*). 

2=1 

The first inequality follows, since E{Tf^p) U 9*}) covers a Steiner tree on {qi,..., qk}. □ 

Figure [2] shows that the factor (2 -|- e) is sharp. For the instance Pi = {(0,0), (1,0)}, P 2 = 
{(0,0), (—1,0)} a minimum two-level rectilinear Steiner tree of length 2 is shown on the left with 
l{Ttop) = 0. On the right, a bad choice of connection points and minimum Steiner trees Ttop,Ti, 
and Pa yield a total length of 4. 
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Figure 2: A tight example when choosing connection points as arbitary points of Pi. 
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Figure 3: A flat Steiner Tree in a lifted instance. 


3 PTAS for a bounded number of partitions 

We can reduce the two-level rectilinear Steiner tree problem in the plane to an ordinary rectilinear 
Steiner tree problem in a higher dimensional space, where we can apply Arora’s PTAS [2]. 

The idea of the PTAS is to lift every subset Pi,..., Pfc to an additional dimension. We assume 
k > 1. Otherwise the two-level Steiner tree problem is an ordinary Steiner tree problem. Let 
Pi,..., Pfc C be the subsets of a two-level Steiner tree instance, we dehne a Steiner tree instance 
in The set of terminals P' is comprised as follows. 

For each original terminal x € Pj C (i € {1,..., k}), we add a terminal := {x, K-ei) G 
where G R^ is the unit vector with value one at the i-ih. coordinate and iL is a large constant, e. 
g. we could choose K as 1{B{P)). Now for x & Ph and y € Pi the distance of their high dimensional 
copies x',y' G P' is ||x' - y'\\i = ||x - y||i 2K\\eh - ei||i = ||x - y||i -b 2K5h,i, where is one if 
h = i and zero otherwise. An example of a lifted two-level Steiner tree is given in Figured 

A (k + 2)-dimensional Steiner tree is called flat if all Steiner points have either the form (x, 0) G 
R 2 -i-fc Qj. (^x,K ■ Ci) G where x G R^ and Cj G R^ is a unit vector. The following Lemma has 

essentially been proven by Snyder m, who shows that an optimum Steiner tree can be found in 
the d-dimensional Hanan grid [9]. We give a short constructive proof for our case. 

Lemma 3.1. A {k + 2)-dimensional Steiner Tree T for P' of length 1{T) can be transformed in 
strongly polynomial time into a {k -\- 2)-dimensional flat Steiner tree T' of length at most l(T). 

Proof. We assume that all vertices and segments of T are located within the bounding box B{P') = 
B{P) X [0, AT]^ of its terminals. Otherwise we could project T into the box without increasing 
the length. Furthermore we assume that all edges in E{T) are one-directional segments after 
introducing vertices of degree two if necessary. 

Fix a dimension j G {3, ... ,k 2} and call an edge in direction j straight. Let X be the set of 
straight segments and let P be the forest that arises from T after removing all straight segments. 
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We construct an undirected graph G as follows. The vertices of G consist of the trees in F plus 
two extra vertices s and t. We add an edge between s and each vertex whose underlying tree is 
located in the hyper plane Hq := {x : Xj = 0}. Similarly we add an edge between the vertices 
whose tree is located in Hj := {x Xj = K} and t. Finally, for each straight edge in X we insert 
an edge between the two vertices representing the corresponding trees in F. 

We assign all edges in E{G) a unit capacity and compute a minimum s-t-cut d{S), where S C 
V{G) with s (z S. is assembled from T by projecting all subtrees in F whose corresponding 
vertex is in S into Hq and all other trees into Hj. To get a valid Steiner tree, we eliminate potentially 
arising cycles by removing edges arbitrarily. By the projection, straight edges in the cut 5{S) will 
grow to a length of K and straight edges outside will shorten to length zero and can be eliminated. 

By Menger’s theorem m there are |<5(5)| edge-disjoint paths between s and t. Each path 
represents a set straight edges of total length at least K. Thus 1{X) > |(j(5)| • K and 

1{T^) = |(5(5)| • K + 1{F) < 1{X) + 1{F) = 1{T). 

Frist, we apply this transformation to all dimensions j G {3,... A: -|- 2} to obtain a Steiner tree 
with Steiner points x of value either xj = 0 or xj = K in the dimensions j G {3,... k + 2} without 
increasing the length of T. Second, we remap all Steiner points x with at least two coordinates 
j, i G {3,... , k + 2} of value Xi = Xj = K to x {0} by setting all Xj j G {3,k + 2} to zero, 
we obtain a flat Steiner tree T' of lenght at most l(T), since we are in the rectilinear case. 

The running time is dominated by the k minimum cut computations that can be done in a 

digraph, where each edge in E{G) is represented by two oppositely directed edges, in strongly 

polynomial time [8]. □ 

Next we show that we can assume that a tree T for P' has at most one edge in direction 
j G {3,..., k}. 

Lemma 3.2. If K > 1{B{P)), a {k + 2)-dimensional flat Steiner tree T for P' of length 1{T) can 
be transformed in strongly polynomial time into a {k + 2)-dimensional flat Steiner tree T' for P' of 
length at most 1{T) so that T' contains at most one edge in each lifting direction j G {3,... A: -|- 2}. 

Proof. Assume that there is a direction j G {3,... , A: -|- 2}, s.t. T contains two edges e, e' in 

direction j. Removing e' splits T into two components which we can re-connect by an edge in 
either Hq := {x : Xj = 0} or Hj := {x : Xj = K} between an endpoint of e and e'. The new edge 
has length at most 1{B{P)) < K, the length of e'. □ 

The proof also shows that if K > l(B{P)) and T has minimum length, it contains at most one 
edge in each lifted direction j G {3 ,.. .k-\- 2}. The following lemma shows the equivalence between 
the original two-level rectilinear Steiner tree problem in the plane and the lifted regular rectilinear 
Steiner tree problem. 

Lemma 3.3. If k > 1, a two-level rectilinear Steiner tree T for Pi,... ,P}^ of length 1{T) can be 
transformed into a {k-\-2) dimensional Steiner tree T' for P' of length at most 1{T) kK and vice 
versa. 

Proof. W.l.o.g. we assume that the center of the coordinate system coincides with the bounding 
box center of P.Let T = {Ttop,Ti,... ,T}f) be a two-level Steiner tree of length 1{T). We embed 
the vertices x G V(Ttop) C into x {0}^ as (x,0), where 0 is a k-dimensional zero vector. For 


5 


Ti (i = {1,..., A:}, we embed x e V{Ti) C as (x, K ■ Ci) ^ R2+^, where e* € is again a unit 
vector with value one in its i-th coordinate. 

To connect the top-level and its subtrees, we pick for each subtree i G {1,..., A;} a connection 
point qi G Ttop H Tj and connect the lifted components by a new edge {{qi, 0), {qi,K ■ Cj)} of length 
K. Clearly, the length of the lifted tree is 1{T) + kK. 

Now let T' be a rectilinear Steiner tree for P' of length 1{T'). By applying Lemma l3.II we obtain 
a (A: -|- 2)-dimensional flat Steiner tree. Applying Lemma 13.21 we can further assume that for each 
j G {3,..., A: -|- 2}, T' contains at most one edge in direction j. 

Removing all k edges of length K, T' is decomposed into A; -|- 1 subtrees. Projecting these onto 
the first two coordinates we obtain a feasible two-level Steiner tree of length at most 1{T') — kK. □ 

Theorem 3.1. For bounded k there is a PTAS for the two-level rectilinear Steiner tree problem. 

Proof. Choose K = 1{B{P)) and for e > 0 set e' := Then compute an (1 -|- e')-approximate 

{k -\- 2)-dimensional Steiner tree T' for the lifted terminal set P' with Aroras PTAS [2] that has 
a polynomial running in bounded dimension. Then we apply Lemma 13.31 to obtain a two-level 
Steiner tree T = {Ttop,Ti,... ,Tfc) for Pi ,... ,Pk with length at most 1{T') — kK. Let T'* and T* 
be optimum Steiner trees for P' and P. Since 1{T*) > K the length of T is 

1{T) < 1{T') - kK < {1 + e')l{T'*) - kK < {1 + e')l{T*) + {1 + e')kK - kK 

< (1 + (A; + l)e')/(r*) = (1 + e)l{T*). 


□ 

Note that [|] +2 dimensions would also be sufficient to achieve this result by lifting two partitions 
into one extra dimension (one partition in positive direction and the other in negative direction). 
This lifting method does not work using Steiner tree approximation algorithms with a constant 
factor a, because we get an additional additive error of a • A: • iL, which in total is essentially not 
better than the 2a-factor approximation by the simple algorithm in Section [2l 

Our construction transforms an optimum (A; -|- 2)-dimensional Steiner tree for P' in the {k 2)- 
dimensional Hanan grid, which always exist due to Snyder m, into a two-level rectilinear Steiner 
tree in the 2-dimensional Hanan grid defined by P. 

Corollary 3.1. Given a set of terminals P = (Pi,... ,Pfc), there is always an optimum two-level 
rectilinear Steiner tree for P in its Hanan grid. 

This also proves the following simple fact. 

Corollary 3.2. The two-level rectilinear Steiner tree problem is in NP. 


4 Predetermined Connection Points 

In all algorithms of this section we predetermine a connection point q^ for each set Pi {i = 1,... ,k) 
and then call a Steiner tree approximation algorithm for {qi,... ,qk} to get Ttop and Pi U {qt} 
to get Ti {i = l...,A:). We use the fact that we consider rectilinear instances to obtain better 
approximation factors than in Section [2l 
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Figure 4: A tight example when choosing connection points as bounding box centers. 

4.1 Bounding Box Center 

A natural approach is to choose each connection point qi (i € as the center of the 

bounding box B{Pi). 

Theorem 4.1. Using bounding box centers as connection points, we get a 1.7f)a-factor approxima¬ 
tion algorithm for the two-level rectilinear Steiner tree problem, when using an a-factor approxima¬ 
tion algorithm for rectilinear Steiner trees as a subroutine. 

Proof. Let T be the resulting two level Steiner tree and let T* = {Tf^p, Tf,... ,T^} he a minimum 
two-level Steiner tree. We choose under all minimum two-level Steiner trees as long as possible 
so that we can choose connection points qf € Tf^p n T* n B{Pi) for all i € {1,..., k}. 

Now, consider a partition set Pi {i G {1,..., k}). We assume w.l.o.g. that the horizontal length 
of B{Pi) is no less than the vertical length. Let q[ be a point in the intersection of T* and the 
vertical line through qi. Note that \\q[ — qi\\i < jl{B{Pi)) by the shape of B{Pi). 

Now for each i e {1,... ,k}, E (T^p) U covers a Steiner tree for {gi,..., and 

E{T*) U {qi,q'i} covers a Steiner tree for Pi U {qi}. Thus, 

k 

1{T) = l{Ttop) + ^l{Ti) 

i=l 

k k 

< a ■ l{Tt,p) + a - qiWi) + l{Ti) 

i=l i=l 

k k k 

< a ■ l{Tfgp) + a ^(||<?* — gdli) + -|- a^(||g' — gi||i) 

2=1 2=1 2=1 

<a-l{Tf,p) + '^-aY,m) < +< 1.75a-/(r*). 

2=1 2=1 


□ 

Figure 0] shows that the factor (7 -|-e) is sharp. For the instance Pi = {(0,0), (1,0), (0,1)}, 
P 2 = {(0,0), (—1,0), (0, —1)} a minimum two-level rectilinear Steiner tree of length 4 is shown on 
the left with l{Ttop) = 0. On the right, with bounding box centers as connection points and shortest 
Steiner trees the total length is 7. 
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4.2 Adjusted Bounding Box Center 

We can improve the approximation factor by a more careful choice of the connection point. For a 
set Pi {i G {1,... , k}), we call the coordinate system with origin in the central point of its bounding 
box the coordinate system of Pi. If no terminal in a subtree instance Pi is located in the lower 
left quadrant of the bounding box w.r.t. its coordinate system, it appears reasonable to shift the 
connection point to the upper right towards the actual terminals, e. g. one would move q 2 towards 
the upper right in Figure 01 

Otherwise, if a set Pi of terminals contains an element in each quadrant of its bounding box B{Pi), 
we call the bounding box B[Pi) complete, as the left example in Figure [5] shows. For subtrees with 
a complete bounding box we choose the connection point to the top-level tree as the central point 
of the bounding box as in Section 14.11 

For sets Pi with an incomplete bounding box B[Pi) we assume w.l.o.g. (after reflection and 
rotation) that there is no point of Pi in the lower-left quadrant of the coordinate system for Pi and 
that the horizontal length is not less than the vertical length of the bounding box B{Pi). We then 
define 


= max{s G : P n {(x,y) : x < s,y < s} = 0}, 

ff = min{s G [0, : P n {(x, y) : x > s, y > s} = 0}, and (2) 

ti = m.\ii{t\,ti,ptf^^}, 

where /3 G [0,1] is a parameter we will choose later and is half the vertical length of B{Pi), 
which by our assumption is also the (vertical) distance from the center to the boundary of P(Pj). 

For subtrees with incomplete bounding box we choose the connection point Qi = {qi^, qiy) to the 
top-level tree as the central point of the bounding box plus the vector {ti,ti) (see also Figure [5]). 
As before, we compute Steiner trees Ttop for {yi,..., qk} and Tj for Pi U {qi] (z G {1,... , k}). 

Moreover for i G {1,... ,k} the trees Ti are refined as follows. We compute a second Steiner 
tree T[ for Pj U {yj} as follows: We compute a Steiner tree T/ for Pi. Thereby, we embed maximal 
paths in containing only Steiner vertices with degree two so that each such path has minimum 
distance to qi while preserving its length. We then add an edge from qi to o^, where o* is a point 
in T] minimizing the distance to qi. Finally, if l{Tf) < l{Ti), we replace Ti by T{. 

Let T be the two-level Steiner tree computed as described and let T* = Tf,... ,T^} he an 




Figure 5: On the left an example for a complete bounding box with q^ chosen at the center. On 
the right an example for an incomplete bounding box. {tj,tj) is the bend of the red line, 
here dehning qi, and (t‘f,t‘f) the bend on the orange line. 

















Figure 6: An example of the situation in the proof of Lemma [4. II The green diamond is the /i-circle 
with radius h around qi. 

minimum two-level Steiner tree. Again we choose under all minimum two-level Steiner trees 
as large as possible so that there is a connection point q* = {qi^,qiy) G Tjop H T* C B{Pi). 

We start with two lemmas bounding the distance \\qi — Oj||i between the connection point qi and 
a nearest neighbor Uj in T/ by the length /(T*) of the subtree in an optimum solution. 

Lemma 4.1. Let i G A;} and T-,ai he constructed as above. If B {Pi) is complete, then 

l{B{Pi)) + \\qi-aih<l{T*). 

Proof. Define h := \\qi — ai||i. Since B{Pi) is complete, T/ intersects at least three of the four axes 
of the coordinate system to Pi. We assume w.l.o.g. that T/ intersects the left, upper and right axis. 
(We did not rotate complete boxes before but only incomplete ones.) 

By the choice of T/ and h there exist (see also Figure [6]) 

p € {(x, y) € Pi : X < -h,y < 0}, 
p' € {(x, y)ePi:x>h,y<0}, 
u € {(x, y) ^ Pi ■. X <D,y >h}, 

Vz '■= {{x,y) ^ Pi '. X > z,y > h — z} for all z G [0, h]. 

Let Vh ^ Vh. If the unique T*-paths from p to u and from p' to Vh intersect, then T* connects 
the lines {(x,y) : x = 0} and {{x,y) : x = h} twice, and therefore l{B{Pi)) + h < l{Tf). 

Otherwise, we can choose a minimum z >0 such that there is a u G 14 and the unique T*-paths 
from p to u and from p' to v are disjoint. The lines {(x,y) : y = 0} and {(x,y) : y = h — z} are 
connected twice in T*. Therefore we get l{B{Pi)) + h — z < 1{T*). 

If z = 0 we are done. Otherwise, if our statement is false there is an 0 < e < z such that 
e = l{B{Pi)) + h — l{Tf) and a w € • Since the unique T*-paths from p to u and from p' to w 

are not disjoint, T* connects the lines {(x, ?/) : x = 0} and {(x, y) : x = z — ^} twice. Therefore we 
get the contradiction 

e = l{B{Pi)) + h- 1{T*) < l{B{Pi)) + h- l{B{Pi))-h + z-z+^-=^-. 

□ 

Lemma 4.2. Let i G k] and T[,ai he constructed as above. If B {Pi) is incomplete and 

ti < then l{B{Pi)) + \\qi - aj||i < l{Tf). 
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Proof. If tj < tj all {x, y) € Pi are either x < t‘f,y > tj or x > tj,y < tf. By the construction of T/ 
we have a* = qi and we are done. Otherwise, there is a point p = {px,Py) G Pi with either px = tj 
and py < tj or px < tj and py = tj. W.l.o.g. px < tj and py = tj. Let be p' € P* a point in the 
lower border of the bounding box. The TLpath from p to p' intersects the left, upper and right 
halfline starting in q^. We are now in the same setting as in the proof of Lemma l4.ll We obtain 
analougsly l{B{Pi)) + \\qi - aj||i < l{Tf). 

□ 


Lemma 4.3. If B{Pi) is incomplete, q*^ < qix, and q*y < qiy, then 

M - 9.111 < jm) + ( 5/3 - 1 ) c“. 

Proof. Since B{Pi) is incomplete, qf € B{Pi), and by the choice of ti and qi in (l2|), there are points 

u G {{x, y) € Pi-. X < q*,^, y > qiy}, 

V G {ix,y) e Pi:x> qix,y < q^y}, 
w G {{x,y) ePi:x> qix,y > qiy}. 

By the positions of the four points qf, u, v, and w, T* connects either the lines {{x,y) : x = qix} 
and {{x,y) G : x = q*^} or the lines {{x,y) G : ?/ = qiy} and {{x,y) : y = q*y} twice. 
Therefore we have minUgj,,, - q*J, \qiy - q*y\} + l{B{Pi)) < 1{T*) and 

Wdi II1 \^i X I T ll/i Qiy I 

< I (.KTf) - l{B{Pi))) + ^mm{\qi^ - q*J, \qiy - qfy\} + ma^{\qi^ - qfj, \qiy - q*y\} 

< - liipm + liti+ c “)+U + \ihoriBiP)) 

< + ^ti - If”- < ^liTf) + Q/3 - tr^, 

where lhor{B{Pi)) is the horizontal length of B{Pi) and we first use < lhor{B{Pi)) to dissolve 

minimum and maximum, and then l{B{Pi)) = + lhcyr{B{Pi)) and ti < □ 


Lemma 4.4. If B{Pi) is incomplete and if q* G {{x,y) G B{Pi) : x > qix or y > qiy}, then 
\\qt-qi\\l<^^liB{Pi)). 

Proof. Since qi G {{x,y) G B{Pi} : x = y} we have ||g* - qi\\i < ^l{B{Pi)). □ 

We are now able to bound the length of each subtree Ti by the length of its corresponding tree 
T* in the optimum solution. 

Lemma 4.5. Let 1 < a < 1.5 he an approximation factor for the rectilinear Steiner tree problem. 
We define 

['ll 111 3 3 1 ll 

f(a) = min max < — a H—, — a H— aB, -a - B H— > . (3) 

^ ^ /3e[o,i] \ 8 4’ 8 8 2 4^ 4/ ^ ^ 

Then, for all i = 1,... ,k we have 


a\\qj -qi\\i+l{Ti)<fia)-l{Tf). 
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Proof. We perform a case distinction corresponding to the Lemmas 14.1114.31 and 14.41 
First case: B{Pi) is complete (as in Lemma l4.ip . 

a\\qt - qi\\i + l{Ti) < ■ l{B{Pi)) + l{Ti) < ■ 1{T*) + {I - ]^a)\\qi - ai\\i 


1 


1 


< -a-l{T:) + -{l--a)l{B{P,)) < _«+- /(!;*) 
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1 


where the second inequality follows from Lemma 14.11 and the third inequality from the fact that 
the Steiner tree T[ insects the vertical line trough the origin inside the bounding box. 

Second case: B{Pi) is incomplete, qf^ < qi^, and q*y < qiy (as in Lemma 031) • 

By the maximality of qf has at least two incident edges in T*. Since there in no point in Pj 
which is lower-left of qi, we can assume w.l.o.g., that q* has exactly one incident edge in the upper 
and one in the right direction. The first one intersects {{x,y) : y = qiy} in {q*^, qiy) and the second 
one intersects {{x,y) : x = qi^} in iqi^,q*y). 

Now we reroute in T* the path starting in {qf^,qiy) via q* to {qix^qty) by a path via qi and 
get a Steiner tree T'* on Pi with liT'*) = l(T*). Since Tf* is a Steiner tree on Pj U {(?i} we get 
l{Ti)<a-l{Tn = a-l{T:). 

With Lemma 031 (first inequality), and 4f™^^ < l{B{Pi)) < 1{T*) we get 

aM-qi\\, + l{Ti) < + 

< + •/(!;*) +a Q/?-i)z(P(P,)) 

Third case: B{Pi) is incomplete and qf G {{x,y) G B{Pi) '■ x > qi^ oi y > qiy} (as in Lemma 03)1 . 
If ti = then 

a\\qt - qiWi + l{Ti) < ■ l{B{Pi)) + l{Tl) < ./(P(P,)) + « •/(i;*) + (1 -/3) 

where the first inequality follows by Lemma [4.41 In the second inequality we use that the extra cost 
of connecting qi in T/ is bounded by the distance (1 — between qi and the upper boundary 

of P(P). 

If ti < then 

a\\qt - <h\\i + l{Ti) < la ■/(B(Pi)) + ;(T,) < la ■ ;(7;*) + (1 - la)||§j - aj||i 
< la.i(2r) + l(l-la)i(B(Fi)) < + 1) ((r*). 


The first inequality follows by Lemma [4.4l the second by Lemma [4.2l and the third since H^i —ai||i < 
^max < lj{B{P,)). □ 
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Theorem 4.2. The two-level rectilinear Steiner tree problem can by approximated by a factor of 


,, , . f 11 1 11 3^3 U 1 

t{a) = mm max < — a H —, — a H—op, —a - p A — 

^ /3e[o,i] 1 8 4’ 8 8 2 4^ 4 

using an a-factor approximation algorithm for rectilinear Steiner trees as a subroutine. 

Proof. Let T be the two-level Steiner tree computeted using adjusted center points as connection 
points and let T* = ... ,T^} be a minimum two-level Steiner tree with as large as 

possible so that all q* € B{Pi) for all i = 1,..., fc. Then 


i=l 


liT) = l{Ttop) + Y,KTi) < a-l{TfJ + J2»h*-Qi\\i + J2^iTi) 

i=l i=l 

k 

< a-liTfJ + f{a)Y,liT*) < f{a)-liT*). 


i=l 


The first inequality follows since E {Tf^p) U (U^=i{9i') Qi}) covers a Steiner tree for {gi, ..., q^} and 
the third inequality follows by Lemma 14.51 □ 

We get an approximation ratio of | when all bounding boxes are complete. 

Corollary 4.1. There is a 2.37-factor approximation algorithm with runtime O(nlogn) for the 
two-level rectilinear Steiner tree problem. 

Proof. A minimum rectilinear Steiner tree on I terminals in the plane can be approximated by a 
factor a = 1.5 computing a minimum spanning tree in the Delaunay triangulation in 0{l log/) time 
[12]. Now we apply Theorem 14.21 where /(a) = ^ < 2.37 is determined by /3 = ^. □ 

Corollary 4.2. There is an 1.63-/actor approximation algorithm for the two-level rectilinear Steiner 
tree problem. 

Proof. We approximate the rectilinear Steiner trees in Theorem 14.21 using Arora’s approximation 
scheme [2]. Choosing 0 < e < 0.003, we get the claimed approximation factor, where /(a) is 
determined by /3 := |. □ 


4.3 Small Top-Level Trees 

The adjusted box center algorithm computes the connection points based on each partition individ¬ 
ually ignoring the structure of the top-level tree. If there is small box containing one terminal from 
each set Pi,... ,Pk, we can give a better approximation ratio independent of the structure of the 
partition sets. We define the top-level bounding box Btop{Pi, ■ ■ ■, Pk) as the smallest axis-parallel 
rectangle containing at least one point qi € Pi for each i = 1,... k. It is a simple exercise to see 
that Btop{Pi ,..., Pk) can be computed in O(n^). 

We choose connection points qi € Pi (1 Btop{Pi,.. ■, Pk) and compute the two-level Steiner 
tree {Ttop,Ti,.. .Tk) by an a-factor approximation algorithm for rectilinear Steiner trees. Let 

{Tfgp, Tf,... T^) be an optimum solution. In [6] it was shown that U{k) := ^ | is an upper 


12 



bound for the ratio of the length of a minimum rectilinear Steiner tree on k terminals and their 
bounding box. Then 

k k 

l{T) = l{Tt,p) + J2KTi) < U{k)-l{Bt,p{P^,...,Pk))+aY,KTt)- 

i=l i=l 


Now if l(Btop{Pi, ■ ■ ■, Pk)) is small, the approximation factor for the two-level Steiner tree is essen¬ 
tially dominated by a. 
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